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ABSTRACT 

We incorporate a model for black hole growth during galaxy mergers into the semi-analytical 
galaxy formation model based on ACDM proposed by Baugh et al. . Our black hole model 
has one free parameter, which we set by matching the observed zeropoint of the local correla- 
tion between black hole mass and bulge luminosity. We present predictions for the evolution 
with redshift of the relationships between black hole mass and bulge properties. Our simu- 
lations reproduce the evolution of the optical luminosity function of quasars. We study the 
demographics of the black hole population and address the issue of how black holes acquire 
their mass. We find that the direct accretion of cold gas during starbursts is an important 
growth mechanism for lower mass black holes and at high redshift. On the other hand, the re- 
assembly of pre-existing black hole mass into larger units via merging dominates the growth 
of more massive black holes at low redshift. This prediction could be tested by future gravi- 
tational wave experiments. As redshift decreases, progressively less massive black holes have 
the highest fractional growth rates, in line with recent claims of "downsizing" in quasar activ- 
ity. 

Key words: galaxies: formation — galaxies: bulges — galaxies: starburst — galaxies: nuclei 
— quasars: general 



1 INTRODUCTION 

In the local Universe, luminous, dusty, merger-driven starburst 
activity has long been suspected to have quasar activity associ- 
ated with it dSanders & Mirabel 19961 1. Some authors find that the 
most powerful Seyfert II active galactic nuclei are usually found 
in galaxies which have had a starburst in the past 1-2 Gyr and 
use this observation to argue that the brightest quasars are asso- 
ciated with galaxy mergers, (Kauffmann et al.^OOS I, whilst oth- 
ers claim that the brightest quasars are hosted in elliptical galax- 
ies which are indistinguishable from the general elliptical popula- 
tion ( [Dunlop et al. 2003) . At high redshift, sources detected in the 
submillimeter are thought to be starbursts (Chapman et al. 20041, 
many are associated with galaxy mergers CSwinbank et al. 2004i 
and many show evidence of active nuclei when probed deeply 
in the X-rays, although it appears that the AGN makes a much 
smaller contribution to the powerful submm flux than the starburst 
( [Alexander et al. 2003l l. 

Black holes (BH) display strong correlations with 
the properties of their host galaxy, particularly those 
of the galactic bulge i Kormendy & Richstone 1995[ 
[Magorrian et al. 1998[ INovak et al. 2006l l. Black hole mass 
is observed to scale with the bulge's B-band luminosity 
( [Magorrian et al. 1998[ [Kormendy & Gebhardt 2001 1, K-band 
luminosity dMarconi & Hunt 20031 [Hiiring & Rix 20 04 ), stellar 
mass dMarconi & Hunt 20031 [Hiiring & Rix 2004) and velocity 



dispersion (Terrarese & M erritt 20001 IGebhardt et al. 2000t . We 
refer to these collectively as the 'Mbh — bulge' relations. It 
has long been theorized that galactic bulges form through galaxy 
mergers dToomre & Toomre 1972) , so it is natural to speculate 
that these events also drive the strong correlation between the 
properties of the bulge and the mass of the black hole. 

There is strong evidence for a link between galactic star 
formation and accretion onto central black holes. The evolu- 
tion with redshift of the global star formation rate and the 
luminosity density of optical quasars are strongly correlated 
( [Boyle & Terlevich 1998) ). At low redshift, the ratio of the global 
star formation rate to the global black hole accretion rate for bulge 
dominated galaxies, SFR/ AIbu is ~ 1000, which is remarkably 
similar to the ratio of A/bh /Af bulge dHeckman et al. 2 004). How- 
ever, it is still an open question whether black hole growth is cor- 
related with all star formation equivalently, or whether its strongest 
relationship is with star formation in bursts. 

The physical conditions in mergers and starbursts are 
amenable to fuelling the accretion of material onto a central 
supermassive black hole. Numerical simulations of galaxy 
mergers have shown that the asymmetrical gravitational 
potential present during the merger is responsible for driv- 
ing gas to the centres of the merging galaxies and of their 
remnant in both major mergers (Mihos & Hemquist 1994b] ) 
and minor mergers ( [Mihos & Hemquist 1994a i. The en- 
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hanced supply of gas to the centre of the galaxy leads to 
rapid star formation and is also available to fuel an AGN 
(Norman & Scoville 1988 , Pi Matteo, Sp ringel & Hernquist 2005| 
[Springel, Di Matteo & Hernquist 2005b l. Furthermore, the forma- 
tion of a dense stellar system with a steep _R^^'*-law potential well 
during a gas-rich merger may help to funnel gas to the AGN at the 
very centre. Starbursts appear to be required in the high redshift 
Universe to explain observations of various galaxy populations 
( |Somerville,Trimack & Faber 2001 1 |Baugh et al. 2005| l. The in- 
creased prevalence of starbursts at early epochs may be responsible 
for the accelerated growth of the most massive black holes towards 
high redshift (e.g. Granato et al. 2004, 2006). 

Theoretical calculations of the growth of black holes 
in the Cold Dark Matter (CDM) cosmology in which struc- 
tures grow through gravitational instability have tended to fall 
into one of three classes: (i) calculations based on the rate at 
which dark matter haloes are assembled, either without any 
treatment of galaxy formation (e.g. Efstathiou & Rees 19881 
IHaehnelt & R ees 19931 IHaiman & Loeb 19981 

[Percival & Mil ler 19991 
[Haiman2004; 



[Wyithe & Loeb 2003 



Koushiappas, Bullock & Dekel 2004 
iMahmood, Devriendt & Silk 2005.) 
very simple estimates of the supply of gas ac- 
onto the BH dVolonteri, Haardt & Madau 2003! 



[Yoo & Miral da-Escude 2004 
or with 

creted 

IBromle y, Somerville & Fabian 2004 
(ii) numerical simulations of galaxy 



Islam, Taylor & Silk 2003 [ 
Libeskind et al. 2006b : 



mergers, which use a mixture of smooth particle hydrodynamics 
and simple recipes to follow the fuelling of a supermas- 
sive black hole dCattaneo et al. 2005al Di Matteo et al. 120051 



Hopkins et al. 2006); (iii) semi- 


analytical modelling of the forma- 


tion of ealaxies and black holes dKauffmann & Haehnelt 2000 


Cattaneo 20011 


Enoki, Nagashima & Gouda 2003 


Menci et al. 2004: 'Granato et a 


. 20041 IMonaco & Fontanot 2005 



[Cattaneo et al. 2005bjl. Recently, the semi-analytical approach 
has been extended to produce models in which the evolution 
of galaxies and black holes are coupled, with energy released 
by accretion onto the black hole either truncating ongoing star 
formation or suppressing the rate at which gas can cool in more 
massive haloes C Granato et al. 2004. Monaco & Fonta not 20051 
IBower et al."2006llCroton et al. 2006t . 

In this paper, we incorporate a model for the growth of 
black holes into the Durham semi-analytical galaxy formation code 
GALFORM dCole et al. 20001 ?). Our prescription for growing black 
holes is tied to galaxy mergers and is similar to the first implemen- 
tation of black hole growth in semi-analytical models by Kauff- 
mann & Haehnelt 12000'). Our starting point is the galaxy forma- 
tion model introduced by Baugh et al. (2005). This was the first 
model to match the observed properties of galaxies in both the 
low and high redshift Universe, following the whole of the galaxy 
population and incorporating a self-consistent calculation of the re- 
processing of starlight by dust. In particular, the model reproduces 
the number counts of Lyman break galaxies and sub-mm sources, 
which are both dominated by starbursts. The success of the model 
is primarily due to an increased level of star formation in bursts 
at high redshift compared with previous models, and the adoption 
of a flat initial mass function (IMF) for stars produced in star- 
bursts. The same model also accounts for the metal content of the 
hot gas in clusters and stars in ellipticals ( [Nagashima et al. 2005 a [ 
[Nagashima et al. 2005b[ and the numbers of Lyman-alpha emitters 
(Le DeUiou et al. 2005, 2006). Since the Baugh et al. model has 
been tested extensively, and, in particular, in view of the success of 



this model in reproducing aspects of the galaxy population which 
are associated with starbursts (and hence bulge formation), we have 
chosen to focus on the predictions for black hole growth and quasar 
activity. 

The use of a semi-analytical model allows us to follow a much 
wider population of objects than is accessible by direct numerical 
simulation. This means that we can follow the demographics of 
the black hole population and explore how black holes acquire 
their mass. The latter is of great importance in view of the recent 
observational evidence suggesting that the most massive black 
holes acquired the bulk of their mass at early epochs and that 
it is the lower mass black holes which are being built up most 
rapidly today. This phenomenon has been termed "downsiz- 
ing" dCowie et al. 20031 ISteffen et al.^003l lUeda et al. 20031 
[Bargeretal. 2005[ [Hasinger, Miyaji & Schmidt 2005| ). At first 
sight, downsizing appears to imply that the growth of black hole 
mass is "anti-hierarchical" and thus incompatible with the CDM 
cosmological framework dMarconi et al. 20041 IMerloni 20041 
IShankar et al. 2004b . We will examine here whether or not such 
downsizing is really a problem for hierarchical models of galaxy 
formation. 

The paper is organized as follows. We provide a description 
of the model in ij2l The model contains one free parameter, which 
we set by matching the 2 = Mbh — bulge relations in ^ where 
we also show that our model is consistent with the evolution of the 
quasar luminosity function. In Sj?) we study the growth histories of 
black holes, separating the contributions from black hole mergers 
and direct gas accretion, and show how the relative importance of 
these channels varies with black hole mass. We predict the evo- 
lution of the A/bh — bulge relations and compare this to data in 
SjS] We demonstrate that we are able to produce downsizing in the 
AGN population in Sj6l We summarize our main results, discuss 
their context and outline future improvements to the model in ^\ 



2 METHOD 

In this section we first give a brief overview of our galaxy forma- 
tion model (§2.1, §2. 2), before explaining how the model has been 
extended to follow the formation of black holes (§2.3). We discuss 
the sensitivity of our model predictions to the mass resolution of the 
dark matter merger trees in §2.4. We end this section with a brief 
description of how a quasar luminosity is assigned to an accreting 
black hole, and present some illustrative results for the quasar lu- 
minosity function at selected redshifts (§2.5). 



2.1 The semi-analytical galaxy formation model 

Our starting point is the model for galaxy formation in the CDM 
cosmology described by Baugh et al. ( i2005i ). As we have already 
pointed out in the Introduction, in addition to giving a reasonable 
match to the properties of galaxies in the local Universe, this model 
also reproduces the counts of sub-millimetre sources and the lumi- 
nosity function of Lyman-break galaxies at high redshift. In both 
cases, the model associates these high redshift objects with galax- 
ies which are undergoing merger-driven starbursts. The success of 
the Baugh et al. (2005) model in reproducing observations linked 
with vigorous starbursts and the formation of spheroids is impor- 
tant for the current analysis. Here we will follow the proposal of 
Kauffmann & Haehnelt (2000) and assume that black hole growth 
is driven by galaxy mergers. For a more exhaustive description of 
the physics and methodology behind the semi-analytical model, we 
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refer the reader to Cole et al. ( I2000t and Benson et al. (2003b). A 
gentler introduction to hierarchical galaxy formation may be found 
in Baugh ("20061). 

We will review the aspects of the model which control the out- 
come of galaxy mergers in the next subsection, and will limit our- 
selves here to more general aspects of the cosmological and galaxy 
formation models. We assume a standard ACDM cosmology, with 
a flat geometry, a matter density f^o ~ 0.3, a baryon density Qh ~ 
0.04, a Hubble constant of Ho — 70 kms^^Mpc^^ and a fluctua- 
tion amplitude specified by as ~ 0.9. The break in the local galaxy 
luminosity function is reproduced by invoking a superwind which 
drives cold gas out of galaxies (?; Nagas hima et al. 2005 a^ ; an al- 
ternative physical mechanism to produce this break is AGN feed- 
back in quasi-hydrostatically cooling haloes dBower et al. 20061 
ICroton et al. 2006V Gas cooling is prevented below z = 6 in low 
circular velocity haloes (vc = 60 kms~^), to mimic the impact 
of the presence of a photoionizing background on the intergalactic 
medium dBenson et al. 2002] l. Baugh et al. l |2005l l adopt a timescale 
for quiescent star formation in galactic discs which is independent 
of the dynamical time, which results in gas rich mergers at high red- 
shift (see their fig. 1). They assume that stars which form in bursts 
are produced with a top-heavy initial mass function; this choice 
has no impact on the predictions presented in this paper beyond 
the high fraction of cold gas forming stars that is recycled into the 
IGM. 

The parameters of the Baugh et al. (2005'! galaxy formation 
model are held fixed in this paper; we do not adjust these parame- 
ters in any way when generating predictions for black holes and 
quasars. This is a clear strength of our approach and choice of 
galaxy formation model. Thus, our results are to be viewed as gen- 
uine predictions of the model. The properties of the quasars and 
active nuclei in our model can easily be related to the properties of 
their host galaxies; such comparisons are deferred to future papers. 



2.2 Galaxy mergers 

Mergers between galaxies play an important role in building up 
the mass and determining the morphology of galaxies. When dark 
matter haloes merge in our model, the galaxies they contain are 
ranked in mass. The most massive one is designated as the 'central' 
or 'primary' galaxy in the new dark halo and the remaining galaxies 
become its satellites. The satellites lose any hot gas reservoir that 
they may have had prior to the merger and any subsequent accretion 
of cooling gas is funnelled into the central galaxy. The orbits of the 
satellite galaxies decay through dynamical friction. If the timescale 
for a satellite to sink to the centre of the halo is shorter than the 
lifetime of the halo, then the satellite is merged with the central 
galaxy at the appropriate time (see e.g. lCole et al. 20()0l ). 

The result of a galaxy merger is determined by two princi- 
pal quantities: (i) the ratio of the mass of the accreted satellite to 
the mass of the primary, /merge = Afsmaiier/A/iarger, and (ii) the 
fraction of the mass of the primary disc which is cold gas, /gas ~ 

Afcold.pr imary rimary i 

where A/disc = Mstars + Mcold. If 

the mass ratio exceeds a threshold /eiup, the merger is termed 'vi- 
olent' or 'major'. In this case, all stars present are rearranged into 
a spheroid, with a radius determined by arguments based on the 
conservation of energy and the virial theorem (see lCoIe et al. 20001 
Almeida, Baugh & Lacey 2007). In addition, any cold gas in the 
merging galaxies is assumed to undergo a star formation burst and 
the stars thus produced are added to the new spheroid. In this pa- 
per, we use the parameters set by Baugh et al., who defined major 
mergers by the threshold /eiup = 0.3. 



In cases where /merge < /eiiip, the merger is termed 'mi- 
nor' . In this case, the stars in the accreted satellite are added to the 
spheroid of the primary, leaving intact any stellar disc present in the 
primary. In minor mergers, the fate of the gas in the merging galax- 
ies depends upon the gas fraction in the primary disc and on the 
value of /merge. If the primary disc is gas rich (if /gas > /gas,burst, 
where, following Baugh et al., we take /gas, burst = 0.75), and if 
/merge > /burst (whcrc Baugh ct al. adopted /burst = 0.05), then 
we assume that the perturbation introduced by the merging satel- 
lite is sufficient to drive all the cold gas, from both the primary and 
the satellite, into the spheroid, where it takes part in the burst. Oth- 
erwise, if in a minor merger the gas fraction in the primary disc is 
small, no burst occurs. Alternatively, if the secondary galaxy is very 
much less massive than the primary (i.e. if /merge < /burst) then 
the primary disc remains unchanged, the accreted stars are added 
to the spheroid and there is no burst, irrespective of /gas. The re- 
finements of the Cole et al. ( 2000) model relating to minor mergers 
were described by Baugh et al. ( 12005 1 , who also set the values of 

/burst, /ellip and /gas, burst- 

During a starburst, we assume that all of the cold gas available 
at the start, AIcom, is processed by the burst and suffers one of three 
fates: (i) It is reheated by supernova feedback and returned to the 
hot ISM. (ii) It is ejected from the dark matter halo by a superwind. 
(iii) It forms long-lived stellar remnants. The mass of gas which 
forms long-lived stellar remnants in the burst, AMstars, depends 
on the feedback prescription used, and is calculated as follows: 



AMst 



(1 



burst) H~ /^burst H~ ysw, burst 



Afe, 



(1) 



For completeness, we now define the parameters in this equa- 
tion (for further details see Cole et al. 2000, Granato et al. 2000 and 
Baugh et al. 2005): 

• /9burst = (Vcire,buige/200kms"^)"^, where Kire.buige is 
the effective circular velocity of the bulge. This quantity gives the 
rate at which cold gas is reheated by supernova feedback in units of 
the star formation rate. This reheated gas is returned to the hot gas 
reservoir and is allowed to recool once a new halo forms (i.e. when 
the halo has doubled in mass). 

• /sw gives the rate of ejection of cold gas by superwinds, in 
units of the star formation rate. This gas is ejected from the dark 
halo and is not allowed to recool. This parameter has the following 
dependence on the effective circular velocity of the bulge: 



/sw — /swO 



for Vcire < Vsw 



/sw — /swO X (V^w/V^ire) for Veirc ^ Vs' 



(2) 
(3) 



The superwind feedback model was introduced by Benson et al. 
(2003b). Baugh et al. (2005) set /swO = 2 and Kw = 200km s-\ 

• -Rburst is the fraction of the mass turned into stars which 
we assume is instantaneously recycled into high mass supernovae 
and returned to the cold phase of the ISM. For the flat IMF used in 
bursts, i?burst = 0.41. 

• e fold is the number of e-foldings over which star formation 
(assumed to have an exponentially declining rate) is allowed to take 
place in a burst. We follow Baugh et al. (2005), taking efold = 3. 



2.3 The growth of black holes in galaxy centres 

The observed correlation between the inferred mass of galac- 
tic central black holes and the properties of their host 
spheroids suggests a common origin for these two classes 
of object (e.g. [Magorrian et al. 1998[ IFerrarese & Merritt 20001 
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log.o ( Mb„ / ( h-' ) ) log,„ ( Mb„ / ( h-' Ms ) ) 



Figure 1. The mass function of black holes, computed using different resolutions for the merger trees of dark matter haloes, as indicated by the legend in each 
panel. Our fiducial resolution is shown by the thick solid line. The left hand panel shows the results for 2 = and the right hand panel shows 2 = 6. 



IGebhardt et al. 2000t . We adopt a model of black hole growth sim- 
ilar to that implemented for the first time in a fully-fledged semi- 
analytical galaxy formation by Kauffmann & Haehnelt ( I2000t . 

We assume that any contribution to the black hole mass from 
processes other than galaxy mergers (e.g. the end products of popu- 
lation III stars, primordial black holes or accretion onto black holes 
from galactic discs or from the hot gas within a halo) is negligi- 
ble compared with the change in black hole mass which occurs 
during merger-driven starbursts or through black hole mergers fol- 
lowing a galaxy merger. Hence, the first galaxies to form in our 
model, when gas first cools into galactic disks, do not contain a sig- 
nificant black hole mass. The first important growth of black hole 
mass is assumed to occur during the first merger triggered starburst. 
Our reason for this choice is that the physics of black hole seed- 
ing is very uncertain, and many mechanisms of black hole seeding 
have been proposed, with widely varying associated seed masses 
(see 'Volonteri 2006 and references therein). The largest mass of 
seed black holes suggested by models in the current literature is ~ 
lO'^ - \(fMQ (e.g. lBromm & Loeb 2003 1 1 Volonteri & Rees 20051 
[Lodato & Natarajan 2006) . These models typically apply only to 
metal-free and highly biased regions at high redshift and probably 
therefore only to the seeding of the the most massive black holes 
(they are motivated by the difficulty in producing ~ lQ^h~^ Mq 
black holes at 2 ~ 6 from smaller seeds, since at high redshifts the 
age of the Universe is not long compared to the Salpeter time). 
Other models for seed black holes predict seeds which are less 
massive than lQ^h~^ Mq, which is our black hole mass resolu- 
tion limit. Therefore, it is reasonable for us to neglect the role of 
seed black holes in our calculations. In practice, we assume that if 
a pre-existing seed is indeed required for supermassive black hole 
formation, then the mass is small enough that it only makes a neg- 
ligible contribution to the mass of the final black hole. Futhermore, 
theoretical considerations (e.g. [Begelman 1978', 'Begel man 2002 [ 
|King 2002l > and observations (e.g. Collin & Kawaguchi 2004[l sug^ 
gest that super-Eddington accretion of mass is possible, and we as- 



sume that this occurs during the early stages of black hole growth, 
so that the mass of any seed does not affect the final mass by Ed- 
dington limiting of mass accretion. 

The mass of black holes is assumed to grow during galaxy 
mergers via two channels, accretion of gas during merger-driven 
starbursts and mergers with other black holes. (Note that in the 
recent model by Bower et al. 120061 additional modes of black 
hole growth are considered: accretion during starbursts triggered 
by disc instabilities, and accretion of cooling gas from quasi- 
hydrostatically cooling haloes.) As discussed in §2.2, we allow star- 
bursts, and thus accretion, in both major and minor mergers. (In 
contrast. iKauffmann & Haehnelt 20001 onlv allowed starbursts and 
black hole accretion during major galaxy mergers.) The two chan- 
nels for black hole growth are as follows: 

Firstly, in starbursts triggered by galaxy mergers, we assume 
that a fraction, -Fbh, of the gas mass which is turned into stars is 
accreted onto the black hole: 



AA/bh = Fbh X AAfsi 



(4) 



where AMbh is the mass added to the black hole and AAfstars 
is the mass of stars produced in the burst after taking into ac- 
count feedback processes that may expel gas from the galaxy and 
the recycling of mass from stars (Eq. 1; see ICole et al. 200^ and 
IGranato et al. 2000i . Typically, we use Fbh < 0.03 (this is ex- 
plained in ^STTJ, and so for simplicity we ignore the depletion of 
the cold gas reservoir by black hole growth when we calculate star 
formation. We assume that the growth of black hole mass is not 
limited to the Eddington accretion rate appropriate to our chosen 
radiative efficiency. 

Secondly, if the merging galaxies already host black holes, 
then we assume that these black holes merge when the host 
galaxies merge. In reality, black holes do not merge instanta- 
neously, but gas-dynamical processes are likely to speed-up black 
hole coalescence in gas-rich mergers i jArmitage & Natarajan 2002) 
and circumstantial observational evidence exists to suggest that 
most binary black holes do merge efficiently, even in gas- 
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poor mergers l |Merritt & Milosavljevic 2005[ >. Since we only con- 
sider binary galaxy mergers with instantanous central black hole 
merging, all BH-BH mergers in our model are binary, and 
we ignore slingshot ejection of black holes from the galac- 
tic centre ( [Saslaw, Valtonen & Aarseth 1974[ l. We also ignore the 
recoil of the merger products of unequal mass black holes 
due to the anisotropic emission of gravitational waves, which 
may lead to the resultant black hole being ejected from the 
galaxy nucleus ( IFitchett 19831 |Volonteri, Haardt & Madau 2003 [ 
ILibeskind et al. 2006t . Recent calculations suggest that most re- 
coil velocities are likely to be in the range 10 — lOOkms"^ 
( [Favata, Hughes & Holz 2004) , and thus unimportant except in 
very low mass galaxies. 

Note that we neglect any loss of mass arising from the radi- 
ation of gravitational waves during the merger of two black holes. 
Such radiation could result in the mass of the merger product be- 
ing less than the sum of the masses of the black holes from which 
it formed dYu & Tremaine 2002 >. This effect is very uncertain, but 
is maximal for equal mass black holes, and even then it is likely 
to be small - approximately 3% or less of the initial mass energy 
for equal-mass non-spinning or Kerr black holes (Baker et al. 2002, 
2004). Since most BH-BH mergers in the Universe have unequal 
mass ratios, the cumulative mass loss by gravitational radiation is 
unlikely to be more than the figure of 20% predicted using the most 
extreme models for gravitational wave loss in individual BH-BH 
mergers dMenou & Haiman 2004b . Therefore, we assume that the 
final black hole mass is the sum of the mass accreted plus the mass 
of the two progenitors. 

2.4 Resolution tests 

The black hole mass down to which our predictions for the proper- 
ties of black holes can be trusted depends upon two factors: the 
accuracy of our prescriptions for handling the physical ingredi- 
ents of our galaxy formation model and the mass resolution of the 
dark matter halo merger trees. The semi-analytical galaxy forma- 
tion model gives a reasonable match to the field galaxy luminosity 
function, including its faint end (?). Further tests of the modelling 
of the phenomena operating in low mass systems are deferred to fu- 
ture work. This leaves the mass resolution of the halo merger trees 
as a numerical parameter that directly influences the properties of 
low-mass black holes. 

In this paper, we use dark matter halo merger trees gen- 
erated using the Monte Carlo scheme described by Cole et 
al. J2000b . Merger trees extracted from N-body simulations 
are, in some respects, more accurate (e.g. |Kang et al. 2005 [ 
|Nagashima et al. 2005c| l. However, a major limitation of the trees 
extracted from simulations is their finite mass resolution. Unpub- 
lished work by one of us (CGL) and work in preparation by Helly et 
al. show that the merger trees in the Durham semi-analytic model 
agree well with the merger trees in N-body simulations. Monte- 
Carlo generated trees can have far superior mass resolution, be- 
cause the whole of the computer memory is devoted to one tree at a 
time, rather than to a large ensemble of haloes within a cosmologi- 
cal volume. This also means we are able to extend our merger trees 
to high redshifts (we use Zstart ~ 20). Also, Monte-Carlo trees 
typically have superior time resolution to those taken from N-body 
simulations. On the other hand, Monte Carlo trees tend to become 
less accurate as the time interval over which the trees are grown is 
increased l lSomerville et al. 2000t . 

Putting this caveat aside, we have performed extensive tests 
of the impact of the choice of resolution of the dark matter merger 



trees on our predictions for the mass function of black holes. The 
results of this convergence study are presented in Fig. [T] for z = 
and z — 6. With improved mass resolution in the merger tree, we 
are able to trace more of the gas which cools in low mass haloes be- 
fore reionization. This is the reason for the odd-looking 'bumps' at 
low BH masses in the z = panel. Our fiducial choice of halo mass 
resolution is 5 x I0^h~^ Mq. This is an order of magnitude better 
than the resolution used in our standard galaxy formation calcula- 
tions, and thirty times better than the resolution of the best N-body 
merger trees currently available within a cosmological volume (the 
Millennium simulation of Springel et al. 2005al which can resolve 
haloes of mass 1.72 X lO^^h-'MQ). With our fiducial halo mass 
resolution, our predictions for the mass function of black holes have 
converged for masses of 10^ AIq and above. 



3 DEFINING THE MODEL : COMPARISON WITH 
OBSERVATIONAL DATA 

We first fix the value of the main parameter in our black hole model, 
Fbh, which determines the mass accreted onto the black hole dur- 
ing a starburst (see Section |23t . In §3.1, we set _Fbh by requiring 
that the model should reproduce the local observed relationship be- 
tween black hole mass (A/bh) and the stellar mass of the bulge 
(A/bulge) in which it resides. We also show the model predictions 
for how black hole mass scales with other properties of the bulge. 
Any viable model of black hole growth should also be consistent 
with the observed quasar population. In §3.2, we briefly describe 
how a quasar luminosity can be assigned to accreting black holes, 
and present some illustrative results for the quasar luminosity func- 
tion at selected redshifts. 



3.1 Setting the model parameter: predictions for the present 
day bulge-black hole relation 

The main parameter of our black hole model is the fraction, -Fbh, 
of the mass of stars formed in a starburst which is accreted onto the 
central black hole (after taking into account gas ejected from the 
galaxy by feedback processes and the recycling of mass in super- 
nova explosions and stellar winds, as described in ^2.11 . We fit the 
value of _Fbh by comparing the model predictions to the observed 
correlation between the mass of galactic central black holes and the 
stellar mass of the bulge component, A/bh — A/buigc, as inferred 
by Haring & Rix ( I2004t . Haring & Rix make a dynamical estimate 
of the stellar mass of the bulge. They compile from the literature 
black hole mass estimates made using a variety of techniques (stel- 
lar, gas or maser dynamics). A review of these techniques and their 
uncertainties can be found in Kormendy & Gebhardt (2001). 

We find that a value of Fbh = 0.022 is required for the model 
to match the zeropoint of the observed Mbh — Afbuigo relationship 
(Fig.|2^). It is important to remember that the normalization of this 
relationship is set by the choice of Fbh- However, the slope and 
scatter are genuine model predictions, and as Fig. |2] shows, these 
predictions are in good agreement with the observations. 

Naively one might argue that, since we have assumed that a 
fixed fraction of the mass of stars formed in a burst is added to the 
mass of the black hole, it is hardly surprising that a tight A/bh — 
Afbuigo relationship results. 

We find in the Baugh et al. model that bursts actually play 
a fairly minor role in the formation of bulge stars. The domi- 
nant channel responsible for building up the mass of present day 
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Figure 2. The relation between black hole mass, Mbh, and a selection of properties of the spheroid of the host galaxy. Each panel shows the correlation with 
a different bulge property: (a) the stellar mass of the bulge; (b) the bulge rest-frame K-band magnitude; (c) the rest-frame B-band bulge magnitude; (d) the 
velocity dispersion of the bulge. The model predictions are shown by the line with en'orbars: the line shows the median and the en'orbars the 10-90 percentile 
spread of the distribution. The observational measurements are shown by symbols, with sources indicated in each panel. 



spheroids is the re-assembly of pre-existing stellar fragments dur- 
ing mergers, not the burst accompanying the most recent major 
merger experienced by the galaxy ( |Baugh, Cole & Frenk 1996] l. 
We find that only 15% of the stellar mass in bulges at redshift zero 
was formed in bursts. The other 85% of the stars in bulges was 
originally formed quiescently, in discs, and no black hole accretion 
is associated with the formation of these stars. Thus, the slope and 
scatter of the Mbh — Mbuige relation are non-trivial predictions of 



the model. Essentially, the Mbh — ^^^buigc relation results from the 
evolution in bulge star formation (and in particular the fraction of 
bulge stars which were formed quiescently) in our galaxy forma- 
tion model. This topic is discussed extensively by Croton ( 12006b . 
The scatter is due to the variation in the fraction of the stellar mass 
of a bulge which was formed quiescently in discs, before being re- 
arranged into the spheroid. The slope originates from how this frac- 
tion varies with stellar mass - as bulge mass increases, the fraction 
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Figure 3. The quasar luminosity function at selected redshifts, as indicated in each panel. The model predictions are shown by lines and the data by symbols, 
with the source indicated in each panel. The different line styles con'espond to different assumptions for how the quasar luminosity depends upon the Eddington 
luminosity of the black hole, as indicated by the legend in the bottom right panel. The data are taken from the following papers: Croom = Croom et al. 2004; 
Fan = Fan et al. 2001; Wolf = Wolf et al. 2003. 



of stars which formed in starbursts decreases, so that black hole ac- 
cretion is associated with a lower fraction of the stars in the bulge. 

Further support for both the galaxy formation model and our 
new model for the growth of black holes comes from examining the 
other relationships between black hole mass and observable prop- 
erties of the galactic spheroid, as shown in Fig.|2](b)-(d). 

In Fig.|2](b), we compare our model predictions for black hole 
mass as a function of the K-band magnitude of the bulge with the 
measurements by Marconi & Hunt (20031 Again, the match is very 
good. K-band magnitudes correlate well with stellar mass. In the 
model, the K-band magnitude depends upon the star formation and 
merger history of the galaxy, taking into account all of the progen- 
itors of the galaxy, its dust content and linear size. Observationally, 
this property is completely independent of the bulge stellar mass es- 
timates based on the velocity profile fitting method used by Hiring 
& Rix ( l2004t . 

It is notable that the scatter in both the Mbh — Mbuige and 
Mbh — Mk relations decreases significantly as the bulge mag- 
nitude gets brighter. A number of factors may contribute to this 
result. As shown in ^ less massive black holes vary far more in 
their formation histories than do larger black holes. Therefore, for 
bulges hosting less massive black holes, there is more scatter in the 
time available for stars to form in progenitor discs before starbursts 
and black hole accretion occur. Stars in larger ellipticals and bulges 
tend to be formed earlier. Once stellar populations exceed a certain 



age, scatter in their ages have only a small impact on colour and 
luminosity. 

In Fig.|2](c) we plot black hole mass against the B-band mag- 
nitude of the bulge and compare the model with a compilation of 
data by Kormendy & Gebhardt (20011. The scatter in this relation- 
ship is the greatest of all four variations on the Mbh — bulge re- 
lations shown in Fig.|2] This is due to the sensitivity of the B-band 
magnitude to the details of the recent star formation history of the 
bulge which can vary considerably between galaxies with similar 
mass black holes. 

Finally, in Fig. [2](d), we compare the model prediction for 
the Mbh — o^buigc relation to data from Haring & Rix ( 12004b . 
We calculate the velocity dispersion directly from the circular ve- 
locity of the bulge, assuming (Jbuigc = 1-1 x T4irc,buigc/v'3 
(see I Almeida et al. 2007 1 for an explanation of the pre-factor). The 
full details of the calculation of 14ii-c,buigo are given in Cole et 
al. ( 12000b . We obtain a reasonable match to the data, reproduc- 
ing the tightness of the relationship, except at the largest veloc- 
ity dispersions. For less massive black holes, our model gives 
Mbh oc ffbuigc which compares well with the Tremaine et al. 
( I2002t estimate of the slope of 4.02 ± 0.32. However, for black 
holes more massive than Mbh ~ 10^'^/i~^Mq, the slope is shal- 
lower than observed, closer to Mbh oc o■b^lg^,. Direct accretion of 
cooling gas from a hot reservoir may help to bring the slope of the 
Mbh — f^buigc relation closer to that observed dBower et al. 2006l l. 

The slope of the AfsH — ''"bulge relation should perhaps be re- 
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garded as one of the less robust predictions of the model, because 
of the complexity of calculating (Jbuigc- This quantity depends upon 
the accuracy of the calculation of the radius of the spheroid. Cole 
et al. (2000) introduced a prescription for computing the size of 
merger remnants, by applying the virial theorem and the conser- 
vation of energy to the progenitor galaxies and the remnant. The 
resulting size of the spheroid is adjusted to take into account the 
self-gravity of the disc and bulge and the reaction of the dark mat- 
ter to the presence of condensed baryons. This step is carried out 
using an adiabatic contraction model. The assumptions behind this 
approach are likely to become less valid as the mass of the spheroid 
increases in relation to the mass of the halo. Almeida et al. (2007) 
tested this prescription against the properties of spheroids in the 
SDSS. Whilst the agreement between the observed and predicted 
Faber-Jackson relation (velocity dispersion - luminosity) is encour- 
aging, the predicted slope is somewhat steeper than is observed and 
the brightest galaxies in the model will perhaps have too large a ve- 
locity dispersion. 

The level of agreement with observations that we find between 
different bulge properties and black hole mass is encouraging and 
suggests that, overall, our model of galaxy and black hole formation 
is on a firm footing. 



3.2 Predictions for tlie luminosity function of quasars 

In order to establish further the credentials of our model, we present 
some illustrative predictions for the evolution of the quasar lumi- 
nosity function. Further assumptions and model parameters are re- 
quired to assign a luminosity to the quasar phase which occurs 
when the black hole accretes gas during a galaxy merger. In this 
section, we give a brief outline of our model for calculating the 
luminosity of the quasar and present some results for the quasar 
luminosity function at different redshifts. These predictions are in- 
cluded here for completeness and to allow comparison with previ- 
ous work (e.g. IKauffmann & Haehnelt 2000l l. We will explore the 
form and evolution of the quasar luminosity function in more detail 
in a future paper. 

There are two basic parameters in our model for quasar lumi- 
nosity: the lifetime of the quasar, tq, and the fraction of the ac- 
creted mass-energy that is turned into the bolometric luminosity of 
the quasar, eq . We assume that the gas available for accretion onto 
the black hole in a galaxy merger is accreted at a constant rate, 
M{t), over the quasar lifetime: 

M(t)=AMBHAq for t < tq. (5) 

(Recall that AA/bh is defined by Eq. 1.) We note in passing that if 
we had instead assumed an exponentially decaying mass accretion 
rate, with a timescale given by tq = O.Stbuigo, giving a mass accre- 
tion rate of (i.e. M{t) = AMexp(-t/(0.5tbuige)) ), this would 
lead to very similar results to those we obtain for a constant mass 
accretion rate. 

The quasar lifetime, fq, is assumed to be directly proportional 
to the dynamical time of the bulge, tbuigc- In the simplest case, 
without imposing any further conditions on the luminosity of the 
quasar, this assumption results in a top-hat light curve: 

LQ{t) ^ eQM{t)c for t < tq. (6) 

When computing the luminosity of quasars, the Eddington 
limit may play an important role. A quasar is said to be radiating 
at its Eddington limit when the pressure of the radiation emitted 
following accretion onto the black hole balances the gravitational 



force exerted by the black hole on new material that is being ac- 
creted. The Eddington limit depends upon the mass of the black 
hole. Physical mechanisms have been proposed which permit mass 
to be accreted at rates which exceed the Eddington limit (see, for 
example, [Begelman 1978) . Here, we show the impact of the Ed- 
dington limit on the luminosity of quasars. We consider four differ- 
ent cases: 

Case (1) No Eddington limit is applied to the bolometric lu- 
minosity of the quasar. 

Case (2) The bolometric luminosity is limited by the Edding- 
ton luminosity corresponding to the black hole mass at the end of 
the accretion episode: 

Lci{t) = max(eqM(t)c', LEdd(Mfln,i) ) (7) 

Case (3) The bolometric luminosity is limited by the Edding- 
ton luminosity corresponding to the black hole mass at the start of 
the accretion episode: 

Lq(t) = max(eqM(t)c^ LEdd(M,tart) ) (8) 

Case (4) The bolometric luminosity is limited by the Edding- 
ton luminosity corresponding to the black hole mass calculated dur- 
ing the accretion episode: 

Lq(t) = max(eqM(t)c^ LEdd(M(t)) ) (9) 

Case (4) is the most realistic estimate of the luminosity after 
applying Eddington limiting. However, there is some uncertainty 
in the evolution of the Eddington limit during the burst, as we do 
not know in detail how the mass of the black hole changes from its 
initial value to the final value. Case (2) corresponds to the maxi- 
mum possible Eddington limit during the accretion episode, being 
set by the final black hole mass. Case (3) is the minimum possi- 
ble Eddington limit during the accretion episode, corresponding to 
the initial black hole mass. We find little variation in the quasar 
luminosity function between these three cases, suggesting that the 
precise growth of the black hole over the accretion episode is unim- 
portant. 

We assume that all visible quasars have identical, flat spectra 
over the range of wavelengths of interest, and that a fraction, /b,, 
of the bolometric luminosity is emitted in the B and bj -bands. We 
adopt /b, = 1/15, which agrees well with Elvis et al. (1994). 
Taking into account the bj-band filter profile, we can calculate a 
magnitude for a quasar from its bolometric luminosity: 

Mb,.vcga = 13.2 - 2.5 X logio(/b., X Lboi/10'"'erg s"'). (10) 

Finally, we assume that only a fraction /visible of quasars 
are detected in optical surveys; the remainder are obscured in 
the optical. We first set /visible = 0.25, i.e. only a quar- 
ter of quasars are visible in optical surveys. This is roughly 
in line with the results of X-ray surveys, which typically find 
that 20-30% of quasars are unobscured in soft X-rays (and 
we presume that this is the same fraction visible in the op- 
tical visible in optical) dUeda et al. 20031 ILa Franca et al. 20051 
IShinozaki et al. 20031 [GilU, Comastri & Hasinger 2007) , although 
this fraction is still somewhat uncertain, and furthermore is likely 
to vary with intrinsic quasar luminosity. We then set the parameters 
tq and eg in order to produce a reasonable match to the Croom 
et al. (2004) measurement of the bj-band luminosity function over 
the redshift interval 0.5 < 2; < 2, as shown in Fig. [3] At higher 
redshifts, we show a comparison between our predicted luminosity 
functions and the 1450 A rest frame quasar luminosity functions 
from the SDSS survey jFan et al. 200 It and the combo- 17 survey 
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dWolf et al. 2003t , using the corrections given in the respective pa- 
pers to convert to the B-band, and applying a further minor correc- 
tion to the 6j band. 

To achieve the best fit, we require that the quasar lifetime, tQ , 
be related to the bulge dynamical time, fbuige, by tq = l.Stbuige 
and that the fraction of accreted mass-energy produced as bolomet- 
ric luminosity be eg — 0.06. A typical bulge has a dynamical time 
of tbuigc of 2 X lO^yr at z = 1, 8 x lO^'yr at z = 3 and 2.5 x lO'^yr 
at 2; = 6, although within each redshift bin, the distribution in 
tbuigo is very broad. As noted by Kauffmann & Haehnelt ( 12000b . 
this redshift evolution in tbuigc helps to reproduce the evolution 
in the quasar luminosity function. Our timescales agree fairly well 
with the Martini & Weinberg (2001) estimate of = 4 x lO^yr at 
2 = 2, which is a typical value. We also note that our adopted ra- 
diative efficiency of 0.06 is consistent with standard disc accretion, 
which is likely to be required for optically bright emission. 

Our simple model does a reasonable job of reproducing the 
observed quasar luminosity function at 2 < 2, but over-predicts 
the luminosity function at higher redshifts. Our basic prediction for 
the quasar luminosity function (shown by the solid lines in Fig [3} 
shows strong evolution with redshift which cannot be described as 
pure luminosity evolution. If the Eddington limit is taken into ac- 
count, then the form of the model predictions changes, particularly 
at bright luminosities, where the abundance of objects is strongly 
suppressed, with the result that the predictions match the data bet- 
ter. The suppression affects more objects at higher redshifts - the 
gas supply then is greater for any given mass of black hole, and the 
dynamical timescales are shorter, leading to higher rates of sup- 
ply for any given mass of available gas. The predicted luminosity 
functions are relatively insensitive to the precise details of how the 
Eddington limit is allowed to influence the quasar luminosity. 



4 THE GROWTH OF BLACK HOLES THROUGH 
ACCRETION AND MERGERS 

We present the bulk of our results in this section. The section is 
quite long, so we list the contents here to help the reader navigate 
through the various topics. Firstly, in ^4A\ we give illustrative ex- 
amples of how black holes acquire mass, tracing the mass assem- 
bly history of two black holes. We then explore the demographics 
of the black hole population: in ^4.2\ we present the predictions for 
the evolution of the mass function of black holes and in ^4.3\ we 
investigate how the black holes are distributed between dark haloes 
of different mass. The next few sections deal with how black holes 
build up their mass. In i]4.4l we show the distribution of progenitor 
masses of black holes, and, in ^4.5\ we address the issue of whether 
the accretion of gas or mergers is the main mechanism for accumu- 
lating black hole mass. We present results for the formation redshift 
of black holes in i]4.6| and for their merger rates in i]4.7| Finally, in 
i]4.8l we compare the amount of baryons locked up in black holes 
with other phases, such as cold gas and stars. 



4.1 Illustrations of black hole growth 

Before concentrating on statistical descriptions, it is instructive to 
show some illustrations of how individual black holes grow in our 
simulations. These examples serve to provide a qualitative picture 
of the model, and to make clear certain definitions and results on 
black hole formation histories that will be of use later on. Note 
that, although space limitations restrict us to only two examples. 



there is, in fact, a rich diversity in black hole formation histories in 
the model. 

The mass assembly history of two black holes is shown in 
Fig. Eland Fig.|5] Fig.|4]shows the central galaxy in a halo of mass 
2.9 X 10^^h~^MQ and Fig|5]shows the central galaxy in a halo of 
mass 8.2 x 10^^ Mq. The main part of each panel follows the 
mass assembly tree. Various components are plotted, as indicated 
by the key at the top of each plot: black hole mass, bulge stars, 
disc stars and cold gas. The area of the symbols is proportional to 
the mass in a given component, with reference areas/masses pro- 
vided at the top of each plot. Galaxies containing black holes are 
linked by solid lines, while galaxies not containing black holes are 
linked to their descendents with dotted lines. The redshifts plotted 
are the output redshifts of the simulation. The left-right positioning 
in the plot is schematic, and has no relevance to spatial positions of 
galaxies within the dark matter halo; the 'main branch' (i.e. the 
most massive progenitor at each merger) is always the right-most 
branch of the merger tree. 

In some cases, the black hole in the 'main branch', which we 
denote the 'main progenitor', may not actually be the most massive 
of all the progenitor black holes at a given epoch, particularly at 
higher redshifts. We have chosen to avoid jumping from one branch 
of the black hole merger tree to another when following the 'main 
progenitor' backwards in time. Instead, we start from the present 
day black hole, find its most massive progenitor, and then build 
up a continuous branch by tracking the most massive progenitor at 
each of ~ 25 — 30 output redshifts. 

In the side panel of each figure, we plot the cumulative masses 
of the black hole and bulge stars as a function of time. The solid line 
shows the total mass in these components, adding together all of the 
progenitors. The dotted line shows the mass in the branch tracing 
back the most massive progenitor of the present day black hole (the 
'main branch'). 

In general, we find a very wide variety of black hole formation 
histories, and we have chosen the ones we plotted to be illustrative. 
The formation trees of the most massive black holes tend to be 
too large and complicated to plot effectively. Meanwhile, there is 
a high abundance of black hole trees with just one burst in their 
history, which were not very interesting to plot. All merger trees 
are included, however, in the quantitative results we present later. 

Inspection of the mass assembly trees, particularly the one for 
the more massive galaxy, reveals that there can be many branches 
to the black hole merger trees at high redshifts. However, most of 
the black hole mass is contained in one or two main branches, as 
shown by the closeness of the solid and dotted lines in the side 
panels. In the Baugh et al. (2005) model, the quiescent star forma- 
tion timescale is independent of the dynamical time. This results 
in discs which are gas rich at early epochs (blue circles), with sig- 
nificant quantities of stars only forming at relatively recent epochs 
(green circles). At later times, it is also apparent that the ratio of the 
stellar mass of the bulge to the mass of the black hole is increasing. 
We will present predictions for the evolution of the AIbh — bulge 
relations in a later section. 

4.2 Black hole mass function 

The black hole mass function at various redshifts is shown in Fig. 
[6] The high mass end advances to higher masses at lower redshifts. 
This is unsurprising in a hierarchical galaxy formation model, and 
reflects the corresponding evolution of the dark matter halo mass 
function. The predicted evolution in the black hole mass function 
is quite strong. This is in contrast with observational claims that 
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Figure 4. An example of a mass assembly tree of a black hole and its host galaxy. Progenitor galaxies without black holes are connected by dotted lines. The 
trees show the relative amounts of cold gas, disc stars, bulge stars and black hole mass, as indicated by the key. The area of the symbols is proportional to mass. 
The left-hand side panel shows the assembly of black hole mass and bulge star's, adding all progenitors (solid line) and tracing back the main branch, which 
is usually the most massive progenitor (dotted hne). The left-right positioning in the plot is purely schematic and has no relevance to the spatial positions of 
galaxies. The final galaxy is the central galaxy of a halo of mass 2.9 X 10^^h~^ Mq. 
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Figure 5. A second example of a mass assembly tree of a BH and its host. The final galaxy is the central galaxy of a halo of mass 8.2 X lO^^h ^Mq. 
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the abundance of large black holes does not vary with redshift (e.g. 
[McLure & Dunlop 2004^ . Such studies, however, typically include 
only optically selected quasars, and so can only probe accreting 
black holes. We examine the relationship of accreting black holes 
to the general black hole population in a later section on downsizing 
(!j6) and consider the implications of the rare, massive black holes 
at high redshift inferred from the observations of Fan et al. ( 1200 It 
in the discussion (0. 

At redshift zero, the break in the black hole mass func- 
tion occurs around l(fh~^MQ. This corresponds to the scale at 
which there is a transition between accretion-dominated growth 
and merger-dominated growth (as we demonstrate specifically in 
i]4.5l ). In larger galaxies hosting more massive black holes, the cold 
gas has been substantially depleted, so the black hole mass can 
only increase significantly through mergers. Gas depletion and sup- 
pression of further cooling by feedback processes is also the likely 
mechanism by which a break in the luminosity function of galaxies 
is produced (?; |Bower et al. 2006llCroton et al. 20061 1. 

From the observed Mbh — Mk, bulge relation, Mbh = 
l(fh~^MQ corresponds to -MK,vega(buige) — Slogig/i ~ —23.5. 
This is very close to the break in the K-band luminosity function, 
M*K - 5 logio/i = -23.44 ± 0.03 dCole et al. 20011 . Similarly, 
from the observed Mbh - Mbuige relation, Mbh = lO^/i'^M© 
corresponds to Mbuige ~ 3 x 10^^ M©. This is the stellar mass at 
which Kauffmann et al. ( 12003b find a transition in galaxy proper- 
ties. Although black hole mass is related to bulge properties only, 
the identification of the knee in the black hole mass function with 
a transition in the global properties of galaxies is reasonable since 
galaxies brighter and more massive than the transition mass tend to 
be bulge-dominated. The conclusion is that galaxies (particularly 
galactic bulges) and black holes grow together (as demonstrated 
graphically in the side panels of Fig.|4]and Fig.O. 

As time advances, the black hole mass function becomes pro- 
gressively flatter at the low-mass end. To a large extent, this is a 
generic feature of a hierarchical mass assembly model in which 
small objects merge into larger objects (at least when this effect is 
not exceeded by the production of new low mass objects). A further 
contribution to the change in slope comes from less massive black 
holes accreting larger amounts of gas as a fraction of their mass 
than larger ones (i.e. downsizing, ij6}. The combination of these ef- 
fects is greater than the effect of the formation of new, lower mass 
black holes; most black holes are seeded at high redshift, as dis- 
cussed in ij4.6l 



4.3 Black hole demography: the conditional mass function 

Black holes of a given mass form in haloes with a broad range 
of masses. The contribution to the black hole mass function from 
different ranges of dark matter halo mass are shown in Fig. [T] At 
the high mass end of each of these conditional mass functions, 
there is a peak and a cut-off. The peak corresponds to the mass 
of the black hole in the central galaxy, which increases as the mass 
of the central galaxy bulge which is strongly correlated with the 
mass of its host halo. This phenomenon is not restricted to black 
holes. Eke et al. J2004| l find tentative evidence for a similar bump 
in the galaxy luminosity function of groups and clusters which they 
attribute to central galaxies, although this remains controversial 
( [Yang et al. 2005) 1. 

In galaxy formation models, a bump is sometimes present in 
luminosity functions where only galaxies in a limited range of halo 
masses are selected - this is because of the contribution of central 
galaxies dBenson et al. 2003at . This reflects the different physical 
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Figure 6. The evolution of the black hole mass function with redshift, as 
indicated by the key. 



processes relevant to central and satellite galaxies: in the model, 
satellite-satellite mergers are not allowed, while all cooling gas is 
funnelled into the central galaxies. These simple assumptions, com- 
mon in semi-analytic models, have been validated in gasdynamic 
simulations ( [Zheng et al. 2004) . Models with intense star forma- 
tion in bursts, such as the Baugh et al. ( 2005 ) model used here, 
smear out the bumps somewhat (Eke et al. 2004), since the bursts 
introduce additional scatter in the properties of galaxies that form 
in haloes of a given mass. 



4.4 Mass function of progenitor black holes 

We now consider the distribution of the masses of black hole pro- 
genitors at different epochs for present day black holes. Fig. [8] 
shows the distribution of progenitor masses at different redshifts, 
for two ranges of black hole mass measured at z = 0. The left- 
hand panel shows the progenitors of z = black holes with 
masses in the range 10^ — 10^'^ Mq and the right-hand panel 
shows the progenitors of z = black holes with masses in the range 
lO" _ lO^'^/i"^M0. The different line types in the plot show the 
progenitor mass distributions at different redshifts, as indicated by 
the key. The distributions plotted are averaged over large numbers 
of black holes with the appropriate present day mass. In both pan- 
els, the z — distribution is naturally peaked around the present 
day mass of the black hole. 

The evolution of these progenitor mass functions from z = 6 
to z = 1 looks remarkably similar to that of the universal black hole 
mass function, albeit truncated at the final z = black hole mass, 
and with an overall normalization which increases with increasing 
final black hole mass. 

The similarity of the form and evolution of the progenitor 
mass functions with those of the universal mass function is re- 
markable. Only at the lowest progenitor redshift plotted (z = 
1) for black holes with present day masses in the range 10^ — 
lO^ -'^/i~^M0 do we see a significant deviation from the form of 
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Figure 7. The contribution to the black hole mass function from black holes 
in dark matter haloes of various masses, as indicated by the key. We show 
the case for black holes contained within any galaxy (thick lines) and in 
central galaxies only (thin lines). The black hole mass function for each 
halo mass is normalized to be the mass function per dex in dark matter halo 
mass. 



ers exceeds that by accretion around a redshift of 1 .7, and growth by 
accretion almost halts after thisQ By redshift zero, the cumulative 
mass assembled by mergers greatly exceeds that assembled by ac- 
cretion. The declining importance of growth by accretion for black 
holes of mass > lO' h^^ Mq reflects the decline in the amount of 
gas available in mergers as more and more of the gas in collapsed 
haloes is consumed into stars. 

In Fig. [Tol we plot the fraction of the mass of black holes 
which, by 2: = (left) and 2 = 2 (right), has been accumulated 
by mergers or accretion onto the 'main branch'. This is shown as 
a function of black hole mass. Fig. [TO] (left) shows that, by red- 
shift zero, low mass black holes have accumulated nearly all of 
their mass by direct accretion onto a single 'main branch', while 
the most massive black holes accumulate 80-90% of their mass by 
mergers of less massive black holes onto the 'main branch'. The 
transition from accretion-dominated growth to merger-dominated 
growth occurs at a 2: = mass of just over l{!ph^^ Mq. Fig. 1101 
(right) shows that at 2 = 2, all black holes, even those more mas- 
sive than Hfh-^MQ, grow predominantly by accretion, although 
there is a contribution from mergers which increases with black 
hole mass. Comparison of the results in Fig. [TO] for 2 = and 
2 = 2 shows that for any given black hole mass, growth by accre- 
tion is more significant for a black hole at 2 = 2 than at 2 = 0, and 
that this difference is greater for the more massive black holes. This 
is consistent with the idea that the luminous growth (i.e. growth by 
direct accretion of gas) of higher mass black holes switches off to- 
wards lower redshifts (see §6). 



the overall mass function. The progenitor mass function in this 
case is rather flat, with fewer low mass black holes and more high 
mass black holes close in mass to the final black hole. The ampli- 
tude of the progenitor mass functions is substantially larger for the 
black holes with present day mass of 10^ — l(f '^h~^MQ: larger 
black holes have a significantly larger number of progenitor black 
holes. This fits in well with our later result ( i]4.5b that less massive 
black holes grow primarily by accretion onto a single main branch, 
whereas black holes larger than 5 x lO' Mq grow primarily by 
mergers of pre-existing black holes. 



4.5 Black hole growth by mergers and accretion 

We come now to one of the principal results of our paper, the man- 
ner in which black holes acquire their mass. There are two distinct 
modes of mass assembly in our model: "accretion", in which cold 
gas is turned into black hole mass during a starburst, and "merg- 
ers", in which existing black holes merge to build a more mas- 
sive black hole. In the accretion mode, mass is being turned into 
black hole mass for the first time, whereas in the merger mode, 
pre-existing black hole mass is being rearranged or reassembled 
into a more massive black hole. In Fig. |9l we plot the fraction of 
the mass in the 'main branch' which is assembled by mergers or 
gas accretion as a function of redshift. We show results for black 
holes in two mass ranges at redshift zero: 
(left) and 10^ - IO^-^/i-^A/q (right). Fig.|9]shows that at high red- 
shifts, growth by accretion dominates over growth through mergers. 
Mergers become increasingly important as redshift decreases, but 
for the 10^ - lO'-^h'^MQ black holes, the cumulative growth by 
mergers never exceeds the cumulative growth by accretion, even at 
redshift zero. However, for black holes of mass 
at 2 = 0, the cumulative growth of their main progenitors by merg- 



4.6 The redshift of black hole formation 

In Fig. [TT]we show the formation redshifts of black holes, binned 
by 2 = mass. Each of the six panels corresponds to a different 
definition of formation redshift. 'Formation' is defined as the time 
when either the main progenitor (right three panels) or the sum of 
all existing progenitor black holes (left three panels) first exceeds 
a given fraction of the final black hole mass. Where the formation 
redshift is defined as the time when the main progenitor first ex- 
ceeds a given fraction of the final mass, we refer to this as the mass 
assembly redshift, since this is the redshift where the stated frac- 
tion of the mass has been assembled into a single object. Where the 
formation redshift is defined as the time when the sum of all pro- 
genitors first exceeds a given fraction of the final mass, we refer to 
this as the mass transformation redshift. This distinction between 
the mass transformation time and mass assembly time for black 
holes is analagous to that between the star formation time and stel- 
lar mass assembly time for the stars in a galaxy. We consider three 
different mass fraction thresholds to define formation times: O.OI 
(top), 0.5 (middle) and 0.95 (bottom). 

When we consider the assembly of 50% or 95% of the fi- 
nal black hole mass (Fig. [TT]- middle-right & bottom-right), we 
see clear hierarchical behaviour; the more massive black holes at 
redshift zero peak in their formation times at lower redshift than 
the less massive ones. This is evidence for the hierarchical assem- 
bly of black hole mass into a single final object. However, when 
we consider the redshift at which 50% of the final black hole 
mass has accreted onto any black hole in the merger tree (Fig. 
nil - middle-left), we see the opposite trend; in the mass range 



This redshift can vary from black hole to black hole - z = 1.7 is the 
redshift where the median growth by mergers exceeds the median growth 
by accretion. 
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Figure 8. The mass functions of progenitor black holes for present day black holes with mass in the interval lO'^ - W'^H-'^Mq (left panel) and 10^ - 
lO^'^?t~^M0 (right panel). The distribution of progenitor masses is plotted at different redshifts. as indicated by the key. The mass functions are generated 
by considering a large sample of black holes at z = 0, and the normalization is chosen so that each progenitor mass function is the mass function per black 
hole at 2 = 0. 




Figure 9. The cumulative growth with redshift of the black hole mass in the 'main branch' divided into the contribution from mergers and accretion. We 
consider large samples of black holes with 2 = masses in the range lO'^ - W-^h-^MQ (left) and 10^ - IO^'^/i-IMq (right). The dotted and dashed 
lines connect the medians of the distribution, while the 10-90 percentiles of the distribution are shown as errorbars. 
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Figure 10. The cumulative fraction of the mass assembled by mergers and accretion, as a function of final blaclt hole mass at 2 = (left) and 2 = 2 (right). 
The medians are connected by lines, and the 10-90 percentile spread of the distribution is shown as an errorbar for each black hole mass. 



more massive black holes display 
a distribution of formation redshifts which peaks at higher redshift. 
When we consider the redshift at which 95% of the final black hole 
mass is accreted onto any black hole in the merger tree (Fig. 1111 - 
bottom-left), we see further evidence of downsizing. Although the 
formation rate peaks at a similar redshift {z ~ 1.2) for all black 
holes in the mass ranee A/bh = 10^ - 10^-^/i"^Mq, the decline 
in the fraction forming per unit time as redshift approaches zero 
is far steeper for more massive black holes within this mass range. 
The observational evidence for "downsizing" refers to the accretion 
of mass onto a particular progenitor, which is accompanied by the 
release of energy. Hence, it is the latter trend which is relevant - as 
black hole mass increases, the redshift when mass is accreted onto 
any progenitor increases. We return to this point in ^ 

We consider now the early growth of the black holes. Fig.ll II 
(top-left) shows the redshift when the first 1 per cent of the fi- 
nal black hole mass has collapsed into any of the branches of the 
merger tree. There is a clear trend for larger black holes to be 
seeded earlier. This is also a form of downsizing. All of the black 
holes in our largest mass bin (Mbh = 10^ - lO^'^fe-^M©) and 
many of those in the next mass bin, A/bh ~ 10* — lQ^'^h~^ Mq, 
are seeded before reionization occurs in the model at 2 = 6. An- 
other interesting feature of Fig. [TT] is that, for almost any defini- 
tion of formation time, less massive black holes have a much wider 
spread in formation times than more massive black holes. 

There is little difference in the distribution of formation times 
of black holes of mass A/bh = 10^ - lO^-^ft-^Af© regardless of 
whether we use a definition which relates to the 'main branch' or 
to 'all progenitors'. This follows from our earlier result that black 
holes in this mass range formed almost exclusively by accretion 
onto a single object, with little contribution from mergers between 
black holes ( ij4.5l Fig.llOt. The differentiation between the different 
definitions of formation time begins to become apparent for black 



holes of mass A/bh = 10^ — IQ^'^h^^ Mq and is increasingly 
more significant as black hole mass increases further. This relates 
to our earlier result that the contribution to the final black hole mass 
from mergers of pre-formed black holes compared to the contribu- 
tion from direct cold gas accretion onto the main branch increases 
strongly with increasing black hole mass ( ij4.5l Fig.llOt. 



4.7 Black hole merger rates 

We show the merger rate per unit time of black holes as a function 
of redshift in Fig. [12] We show this for a number of mass thresholds 
which must be exceeded by both of the black holes that take place 
in the merger. The merger rates peak at lower redshift for more 
massive black holes, with the merger rate for the most massive bin 
still rising at « = 0. This is consistent with the trend seen at 2 = 
that larger mass black holes grow primarily by mergers, while less 
massive black holes grow primarily by accretion (^43]l. 

This behaviour in the growth and merging of black holes 
of varying mass is largely a reflection of the general hierarchical 
growth of structure, moderated in the case of galaxies and black 
holes by baryonic processes. The results presented in this subsec- 
tion concern only mergers of black holes, not necesarily their total 
growth which can also involve accretion. There is no evidence for 
'anti-hierarchical' behaviour in the evolution of black hole merg- 
ers. However, as we show in §6, this is perfectly compatible with 
quasar downsizing - black hole merging can be a 'dark' process 
in which no gas is present, whereas the observational evidence for 
downsizing refers to processes involving star formation or gas ac- 
cretion. 
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Figure 11. The distribution of formation redshifts of blacic hole.s in 5 different bin.s of 2 = mass, as indicated by the icey. The differing definitions of 
formation redshift used in each of the six plots are noted briefly on each plot and explained more fully in i|4.6l of the text. 
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Figure 12. Black hole merger rate per unit time as a function of redshift. 
The merger rate is plotted for 5 different mass thresholds, as shown by the 
key, which the (pre-starburst) masses of both black holes must exceed. 



Figure 13. The evolution of the fraction of the baryons in the universe in 
hot gas, cold gas, stars (disc plus bulge) and black holes. The small abrupt 
changes in some of the lines are due to simulation runs that finish at different 
times. 



4.8 The fraction of baryons in black holes 



Having considered the formation of individual black holes, we now 
look at the global picture. In Fig.[T3] we show the integrated cos- 
mic density of all the baryonic components of the universe; hot 
gas, cold disc gas, stars and black holes. After 2 ~ 4, the growth of 
black hole mass in the universe slows down in comparison to that 
of stars, as quiescent star formation begins to dominate over star 
formation in bursts. The decline in cold gas from redshift 2 to 
goes a small way towards explaining the decline in quasar activity 
over this redshift interval. The decline in the galaxy merger rate and 
the transition from burst-dominated star formation to quiescent star 
formation also play a role. In Fig. [141 we show the star formation 
rate, divided into burst and quiescent modes, and the rate of black 
hole growth. By construction in our model, black hole growth is 
more strongly correlated with the star formation rate in bursts than 
with star formation in general. Very broadly, although perhaps less 
so at low redshifts, black hole accretion tracks the overall star for- 
mation over cosmic time, as observed ( [Boyle & Terlevich 1998^ . 

The cosmological mass density of black holes dX z = Q 
is a quantity of interest. In our model, we find that pan = 
2.83 X IO^'A/qMpc"^. Observationally, pen is determined by 
integrating the black hole mass function which, in turn, is in- 
ferred from a combination of the velocity dispersion distribution 
of galaxies, the K-band luminosity function or the bulge stel- 
lar mass function, and the appropriate A/bh — bulge relation. 
Observed values of /9BH/(lO^Af0Mpc~''), converted to Ho = 
70kms"^Mpc"\ are : 2.9 ± 0.5 jYu & Tremaine 2002t . 2.4 ± 
0.8 ( lAUer & Richstone 2002b . 2.8±0.4 ( [McLure & Dunlop 2004b , 
4.2 ± 1.1 JShankar et al. 2004b and 4.6^i:4 ("Marconi et al. 2004). 
Our estimate is towards the lower end of the broad range spanned 
by the observational estimates. We do not include any measurement 
errors in our estimate. A detailed comparison would need to take 
into account galaxy type (some estimates are based only on ellipti- 
cals), the flux limits of the observational samples and the treatment 
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Figure 14. The variation with redshift of the global star formation rate 
(starbursts, quiescent and total) and the global rate of black hole growth. 



of the dispersion in the A/bh — bulge relations when converting 
from bulge properties to black hole masses. Many observational es- 
timates assume that the scatter in log(MBH), at a given value of the 
bulge property under consideration, is symmetrical. This assump- 
tion then leads to larger values of pan for larger assumed values of 
the scatter ( [McLure & Dunlop 2004) . However, it is not at all clear 
that the scatter in these relations is symmetrical. 
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Figure 15. The redshift evolution of the relations between central black hole mass, Mbh, and bulge properties. Each panel shows the relationship between 
A/bh and a different property of the host spheroid : (a) the bulge rest-frame B-band magnitude; (b) the bulge rest-frame K-band magnitude; (c) the stellar 
mass of the bulge; (d) the velocity dispersion of the bulge. The model predictions are shown by the symbols with errorbars; the lines show the median relations 
and the errorbars the 10-90 percentile spread of the distributions. Redshifts 0, 1, 2, 3, 4 and 6 are shown in different line types, as indicated by the key. 



5 THE EVOLUTION OF THE RELATION BETWEEN 
BLACK HOLE MASS AND BULGE PROPERTIES. 

In this section, we discuss the evolution of the relationship between 
black hole mass and various galaxy bulge properties: K-band and 
B-band bulge magnitude, bulge stellar mass and bulge velocity dis- 
persion. We show these relationships in Fig.fTS] in each case plot- 
ting the model predictions for the A/bh — bulge relationships at 
2 = 0, 1, 2, 3, 4 and 6. We discuss each of these in turn, briefly 



referring to any relevant observational data. However, it is difficult 
to make rigorous comparisons to current data. While relationships 
between black hole mass and bulge properties are fairly well de- 
termined at 2 = 0, this is currently not the case for z > 0, where 
observational samples are small and subject to selection effects. In 
particular, different surveys sample the population of galaxies and, 
where relevant, the AGN subpopulation, in ways that are not always 
straightforward to replicate in the models. 
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We show our model predictions for the Mbh — Mb. bulge re- 
lation in Fig. llSf a) and the A^bh — A/k, bulge relation in Fig. llSf b). 
In the model, these relationships shift towards brighter magnitudes 
at higher redshifts. This is a reflection of the evolving stellar pop- 
ulations. The stellar populations of bulges at low redshift are older 
and thus less luminous than their high redshift counterparts. This 
effect more than compensates for any evolution in the opposite di- 
rection in the A/bh — A/buige relation, which we discuss below. The 
redshift evolution in the A/bh — A/b, bulge relation is greater than 
that in the A^bh — Afx, bulge relation because stellar populations 
dim more strongly with time in the B-band than in the K-band. Ob- 
servationally, however, Peng et al. (2006), selecting high redshift 
quasars, find little trend in the A^bh — A/r, bulge, rest relation with 
redshift, which conflicts somewhat with our prediction of an evolu- 
tion towards brighter magnitudes as redshift increases. 

We show the A/bh — Afbuige relation in Fig. llSf c). There is 
no significant evolution in either the slope or scatter at large bulge 
masses. For Afbuige < IQ^^hT^ Mq, the black hole mass to bulge 
mass ratio increases with increasing redshift. Observationally, Peng 
et al. ( 12006b find that the ratio of A^BH/A/buige was 3-6 times larger 
at z > 2 for AGNs than for quiescent galaxies at 2: = 0. McLure 
et al. J200 6 1, selecting radio galaxies at 2 > from the 3CRR cat- 
alogue, argue that A/Bn/Afbuige increases with redshift, and is ~ 4 
times greater for radio galaxies — 2 than for quiescent galaxies 
at 2: = 0. We find that A/Bn/A/buige was ~ 2 times greater at 2 = 2 
than at 2 = for Afbuige < 10^'' /i~Mf0. This evolution is in the 
same sense as and of comparable size to the observational trend, al- 
though the effect in the model is perhaps not as strong. As discussed 
in i ]3.1l the predicted variation in the A/Bn/A^buige ratio reflects the 
variation in the fraction of bulge stars which formed quiescently in 
discs. Mergers at higher redshift, when discs are more gas-rich and 
have fewer stars, deposit a lower fraction of (quiescently-formed) 
disc stars in the bulge (e.g. lCroton 2006t . 

Close inspection of Fig. fTsT c) shows that a few objects at the 
highest redshifts have black hole masses that exceed Fbh x A/buige. 
This would appear to be impossible given our definition of Fbh 
in i]2.3l This apparent anomaly is due to our assumption that the 
mass of the black hole increases instantaneously at the time of 
the starburst. In our model, star formation in the bursts extends 
over ~ 50 dynamical times, while quasars shine only over ~ 0.3 
dynamical times, so that the stellar mass builds up much more 
slowly than the black hole mass. It seems likely, however, that a 
black hole will still be growing towards its final mass towards the 
end of the starburst dArchibald et al. 20021 [Alexander et al. 2005'; 
|Borys et"al. 2005 1. We defer a study of the co-evolution of the stars 
and the the black hole mass to future work. 

Finally, we show the Mbh — crbuige relation in Fig. I15f d). 
There is no evolution in the slope of the relation, but the zeropoint 
does evolve and the scatter increases significantly towards higher 
redshift. For a given mass of black hole, the velocity dispersion of 
the bulge is greater at higher redshift. To some extent, this evolu- 
tion reflects the expected variation in the properties of dark mat- 
ter haloes: at a given mass, the halo velocity dispersions scales as 
(T oc (2form + 1)^'''^. Alternatively, the evolution could be viewed 
as a reduction in the black hole mass with increasing redshift, for a 
fixed bulge velocity dispersion. 

Shields et al. J2003I I have compared the relative amounts of 
black hole mass in distant quasars and in galaxies in the local Uni- 
verse. They find a large scatter and an increase of — 0.5 dex in 
A#BH/a"buige between 2 = and 2 ~ 3. Similarly, Woo et al. 
( |2006t have compared Seyferts at 2 — 0.36 with galaxies at 2 = 0. 
They too find an increase, of 0.62 ± 0.10 dex, in black hole mass 



at fixed (Jbuige at 2 = 0.36 compared to 2 = 0. Thus, the ob- 
served trend in AfBH/o"buige, if any, is in the opposite direction to 
the trend we find in our simulations. It is possible that our model 
neglects effects that would cause black holes to be a larger fraction 
of the galactic bulge mass at higher redshifts. However, it must be 
remembered that (Tbuige is one of the more uncertain properties of 
the galaxies in our model and that dynamical effects which are not 
included could play a role in determining the properties of merger 
remnants dPekel & Cox 2006. .Robertson et al. 2006J . 



6 DOWNSIZING IN A HIERARCHICAL UNIVERSE 

In cosmology, "downsizing" is an ill-defined term which has been 
applied to describe the phenomenon whereby luminous activity 
(e.g. star formation or accretion onto black holes) appears to be oc- 
curing predominantly in progressively lower mass objects (galax- 
ies or BHs) as the redshift decreases. Claims of downsizing were 
first made in connection with the population of star-forming galax- 
ies l lCowie et al. I996] >. More recently, the same trend has been 
inferred from the evolution of the X-ray luminosity function of 
quasars dCowie et al. 20031 ISteffen et al.^003l lUeda et al. 20031 
|Barger et al. 2005[[Hasinger, Miyaji & Schmidt 2 005): the number 
of bright X-ray sources peaks at a higher redshift than the num- 
ber of faint X-ray sources. The optical quasar luminosity function 
shows similar evolution, with more bright objects seen at increasing 
redshifts (e.g. Croom et al. 2004 ). 

The apparent downsizing in the quasar X-ray luminos- 
ity function has been interpreted by some authors as imply- 
ing that black holes acquire mass in an 'anti-hierarchical' 
maimer JMarconi et a l. 2004; Merloni 20041 IShankar et al. 20041 
[Hasinger, Miyaji & Schmidt 2005) . In this Section, we demon- 
strate that the "downsizing" of the luminous growth of black 
holes is actually a natural feature of our model, despite the fact 
that the overall assembly of mass into black holes is hierarchi- 
cal. Downsizing in the galaxy population in hierarchical mod- 
els is promoted by the earlier collapse and more active merg- 
ing of objects in regions of high overdensity (Kauffmann 1995) 
|Mouri & Taniguchi 2006l|Neistein, van den Bo"sch & Dekel 2006) . 
In recent models of galaxy formation, this natural trend is accentu- 
ated by the feedback processes associated with AGN activity in 
massive haloes (Bower et al. 2006; Croton et al. 2006). However, 
we wish to emphasize that AGN activity in low redshift cooling 
flows is very far from being the only ingredient required for down- 
sizing, and that we still find downsizing in our model. We now re- 
view some of the indirect evidence already presented in support 
of this conclusion ( i]6.1b , and go on to present explicit predictions 
which reveal which black holes in our model are accreting mass 
most rapidly ( i]6.2t . 

6.1 Indirect evidence for downsizing in the model: the 
evolution of the optical luminosity function 

The optical quasar luminosity function, as we have already re- 
marked, reveals a dramatic increase in the space density of bright 
quasars with increasing redshift dCroom et al. 2004t . In ^3.2\ we 
presented the model predictions for the optical luminosity function, 
which are in good agreement with this trend in the observations. 
Two features of our model are responsible for this success: the in- 
crease in the halo merger rate (and hence the galaxy merger rate) 
with increasing redshift and the increase in the gas content of discs 
with increasing redshift (see also lKauffmann & Haehnelt 2000b . In 
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combination, these phenomena lead to an increase in the frequency 
and strength of starbursts with redshift. In our model, a starburst re- 
sults in a luminous phase of growth of the supermassive black hole; 
a fraction _Fbh of the cold gas which is turned into stars during the 
burst is accreted onto the black hole. 

Galaxy mergers are still an important way of building black 
hole mass at low redshift. Our model predicts that BH-BH merg- 
ers are the most important channel for building black hole mass for 
the most massive black holes at the present day. This dark growth 
process represents the assembly of mass which is already locked 
up in black holes into larger units. Galaxy mergers at low redshift 
tend to be gas poor in our model simply because more time has 
elapsed to allow galactic discs to turn cold gas into stars quies- 
cently. This effect is accentuated in the case of the most massive 
black holes which tend to reside in the more massive dark haloes. 
The process of galaxy formation starts earlier in the progenitors of 
massive haloes, since these objects collapse into bound structures 
earlier than is the case in less extreme environments. 

6.2 Direct evidence for downsizing in the model: which black 
holes are accreting mass? 

Our model allows us to separate the mass assembly of black 
holes into two contributions: accretion, in which cold gas is turned 
into black hole mass in a starburst and mergers, in which exist- 
ing black holes merge to build a more massive black hole. Here 
we focus on the process of gas accretion. Fig. [16] presents two 
views showing which mass of black holes are accreting material 
the most vigorously. The left-hand panels of Fig. [16] show the 
distribution of accretion rates, expressed in units of the Edding- 
ton mass accretion rate. Since the Eddington accretion rate scales 
with mass, this is easily scaled to give the distribution of frac- 
tional accretion rates. The right-hand panels compare the present 
accretion rate to the past average accretion rate (calculated as 
(M) X tagc of Universe/ (A/bh), as a functiou of black hole mass. 
Each row corresponds to a different redshift (top: z = Q, middle: 
z = 1, bottom z — 2). In these plots, we have not limited the mass 
accretion rate to be less than or equal to the Eddington limit. 

The left-hand side of Fig.[T6]shows that, at all redshifts, there 
is a large spread in the Eddington ratios at which black holes are 
accreting. There is variation amongst mergers in gas supply, accre- 
tion timescale and initial black hole mass. Furthermore, the Edding- 
ton ratio evolves during any single accretion event. As expected, 
the mass accretion shifts towards higher fractions of the Eddington 
limit at higher redshift since there is more gas available in mergers. 
At z = 0, we see that as black hole mass increases, accretion shifts 
to lower fractions of the Eddington ratio. This trend is less pro- 
nounced ai z — 1 and practically disappears by z — 2. Thus, more 
massive black holes were accreting mass more rapidly at z — 2 
than they are today. The predicted distribution of mass accretion 
rates at z = agree reasonably well with the observational results 
of Heckman et al. (20041. The distribution is normalized to unit 
area, but most black holes, particularly at lower redshifts, are not 
accreting at all (i.e. they are in a 5-function at M — 0). 

The right-hand column of Fig. [16] shows the ratio of 
the present accretion rate to the past average accretion rate 
( (M) X tage/ (Mbh) ) as a function of black hole mass at 2 = 
0, I and 2. If this ratio exceeds unity, then the current mass accre- 
tion rate exceeds the average rate at which the black hole gained 
mass in the past (summed over all progenitors). The predictions 
for this ratio are sensitive to the black hole selection, for example, 
selection using a cut in quasar luminosity. We show results for all 



black holes (solid lines) and also for black holes selected as quasars 
brighter than a given luminosity. Note that the bulk of black holes 
in the model are not accreting material at any given time. The solid 
lines in each panel show that there is clear evidence for downsizing 
in the model. At 2 = 0, more massive black holes are growing less 
rapidly than less massive black holes. By 2 = 1, this trend is greatly 
diminished, and at 2 = 2 it is reversed, i.e. the most massive black 
holes have the highest fractional accretion rate. When only quasars 
are selected, those with low mass black holes show similar frac- 
tional accretion rates as redshift varies from to 2, while quasars 
with massive black holes show a strong decline in fractional accre- 
tion rate towards the present day. 



7 SUMMARY AND DISCUSSION 

We have described an extension to the GALFORM semi-analytical 
model of galaxy formation in the ACDM cosmology to track the 
growth of black holes (BH). Our model for black hole growth has 
one free parameter, Fbh, the mass accreted onto the black hole as 
a fraction of the stellar mass produced during a starburst. We set 
the value of _Fbh so as to reproduce the zeropoint of the present 
day A'/bh — bulge relations. The slope, scatter and evolution of the 
A^BH — bulge relations are model predictions. 

In our model, black holes grow only during and following a 
galaxy merger. They grow through two distinct channels: mergers 
of pre-existing black holes and accretion of cold gas if a starburst 
is triggered by the merger. The importance of growth through black 
hole mergers increases with the mass of the black hole; at 2 = the 
growth of black holes less massive than 5 x W"^ Mq is domi- 
nated by accretion, while the growth of more massive black holes 
is dominated by mergers. In general, the growth of black hole mass 
by mergers becomes more important at low redshifts as the sup- 
ply of gas available for accretion is consumed by star formation. 
Our model neglects black hole growth from gas accreted directly 
in a cooling flow from a hot gas reservoir, as may be expected in 
massive haloes at late times. This is the "feedback mode" of black 
hole growth invoked by Bower et al. ( 12006b in their model that 
explains why there is an exponential cutoff at the bright end of 
the galaxy luminosity function. Apart from this new growth chan- 
nel and an explicit treatment of disc instabilities, the calculation 
of black hole growth in the Bower et al. model is very similar to 
ours. However, around 20% of the global mass in black holes in the 
Bower et al. study is due to the "feedback mode" of growth. It is 
as important as mass assembly due to galaxy mergers for the most 
massive black holes, which accumulate ~ 50% of their mass via 
this channel (Richard Bower, priv. comm.). 

Essentially all current observational estimates of the accumu- 
lation of black hole mass are sensitive to luminous growth, i.e. 
mass accretion. However, we predict that the importance of growth 
through BH-BH mergers grows with decreasing redshift and with 
increasing black hole mass. BH-BH mergers represent a dark mode 
of growth that is more difficult to observe and confirm. The most 
obvious way to detect BH-BH mergers is through the emission of 
gravitational waves (e.g. IHaehnelt 1994l l. This may be possible in 
ten years with the planned LISA gravitational wave interferometer. 
When gas is present during a black hole merger, a circumbinary 
accretion disc could form and the BH merger may produce high 
velocity gas outflows ( [Armitage & Natarajan 2002) followed by an 
X-ray afterglow which could be detected by the next generation 
of X-ray observatories ( [Milosavlje vic & P hinney 2005 [ l. Winged or 
X-shaped radio sources ( [Merritt & Ekers 2OO2F and cores in ellip- 
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Figure 16. Left panels: The distribution of accretion rate normalized by the Eddington mass accretion rate. Right panels: the current mean mass accretion rate 
normalized by the past average mass growth rate, plotted against black hole mass. In both cases, each row corresponds to a different redshift: z = (top), 
2 = 1 (middle) and 2 = 2 (bottom). In the left hand panels, each line shows the distribution of accretion rates for black in the mass interval shown by the key. 
In the right hand panels, the lines correspond to different cuts on quasar luminosity, again as shown by the key. 
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tical galaxies dFaber et al. 19971 |Milosavljevic et al. 2002| ) may be 
indirect evidence of gas-poor mergers. 

Our model predicts that the most important growth mech- 
anism for the most massive black holes is the "dark" mode or 
mass assembly through BH mergers. A testable prediction of our 
model is that a tail of black holes with masses above a few times 
IQ^Mq should be found at 2 = once high quality observations 
covering a large volume of the local Universe become available 
(ED- Furthermore, we expect that these will be more massive than 
any found in quasars at high redshift. To date, only black holes 
less massive than ~ 3 x 10^ have been unambiguously ob- 
served in galaxies at 2 = dTremaine et al. 2002) and in luminous, 
optically-selected quasars over the redshift interval < 2 < 2 
( [McLure & Dunlop 2004[ l. However, this implied limit on black 
hole mass is far from robust. The most massive black holes at 2 = 
tend to reside in massive and hence rare elliptical galaxies, which 
could easily have been missed in existing surveys. Larger volumes 
(> a few times Hf'h^^Mpc? ) need to be surveyed to find such 
objects, which are therefore likely to lie at large distances. This, 
coupled with their expected low surface brightness (more massive 
ellipticals tend to have lower surface brightness cores), could make 
it difficult to measure their central mass using methods based on 
stellar dynamics ^Kormendy & Gebhardt 200T) . The high redshift 
quasar data do not give a complete census of the black hole popula- 
tions, as quasar observations are only able to probe accreting black 
holes. 

There is an important distinction to be made in our model 
between mass transformation and mass assembly. Mass transfor- 
mation refers to the process of turning cold gas into black hole 
mass; at any one time in the formation history of a black hole, 
this phenomenon could be occuring across a number of progeni- 
tor black holes. Mass assembly refers to the accumulation of mass 
in a black hole's main progenitor, and may occur via both direct 
accretion of gas and merging of pre-existing black holes. Black 
hole mass assembles hierarchically; more massive black holes are 
assembled at lower redshifts than less massive black holes. How- 
ever, if we choose to define the formation time of a black hole 
in terms of the mass transformation redshift when some fraction 
of its mass has been accreted onto any progenitor, we find that, 
for Mbh > lO^/i"^M0, more massive black holes form ear- 
lier. This dichotomy mirrors the growth of stellar mass in galac- 
tic spheroids in hierarchical models. In the semi-analytical models, 
galaxy mergers produce spheroids. At high redshift, the mergers 
tend to be gas rich and new stars are produced as a result of the 
merger event. At low redshift, galactic discs tend to be gas poor 
and consist mainly of stars, with the result that the merger sim- 
ply rearranges the pre-existing stars (Baugh, Cole & Frenk 1996| 
IKauffmann l"996llBell et ai720 06 ; De Lucia et al. 2006t . 

While we find that black hole mass is assembled hierar- 
chically, our model clearly exhibits a "downsizing" in the mass 
of black holes which are undergoing luminous accretion. At the 
present day, we find that low mass black holes are accreting mate- 
rial at a higher proportion of their Eddington luminosity than high 
mass black holes. This distinction is less apparent at higher red- 
shifts. Another way to demonstrate this downsizing is to examine 
the rate at which black holes are accreting mass, expressed as a 
fraction of the mass already in place: {i\/BH)/(MBH). At 2 = 0, 
A/bh/Mbh is largest for low mass black holes and drops rapidly 
with increasing mass. This trend is removed with increasing red- 
shift; accretion becomes an increasingly important mode of mass 
assembly for all masses of black hole at earlier epochs in the uni- 
verse. 



A number of authors have claimed that black holes grow in 
an "anti-hierarchical" fashion dMarconi et al. 20041 IMerloni 20041 
Shankar e t al. 2004l l. This conclusion is reached by comparing an 
inferred present day black hole mass function with the black hole 
mass function expected from AGN relics under the assumption 
that black holes grow exclusively by accretion. These calcula- 
tions ignore any contribution to the mass of black holes arsis- 
ing from BH-BH mergers. Furthermore, the assumption that all 
black holes accrete at a constant fraction of their Eddington ra- 
tio (Marc oni et al. 20041 IShankar et al. 2004t . or the use by Mer- 
loni (2004) of the 'fundamental plane of black hole activity' 
(Merloni, Heinz & Di Matteo 2003) , which has a very large scat- 
ter, will introduce errors that may become cumulatively very large 
as the black hole mass function is integrated backwards in time. 

While the zeropoints of the present day Mbh — bulge rela- 
tions are set by adjusting a single free parameter, the slope, scat- 
ter and evolution of these relations are genuine predictions of the 
model. We find little evolution with redshift in the slope of any of 
the A/bh — bulge relations, although our model predicts differ- 
ing evolution in their zeropoints, depending upon which particular 
bulge property is being considered. If we focus attention on a fixed 
black hole mass, we find that with increasing redshift, the typical 
host bulge is more luminous in the rest-frame B-band and K-band, 
shows little change in stellar mass (except in the case of low mass 
black holes, where the stellar mass is lower) and has a somewhat 
higher velocity dispersion. 

Our model predicts the presence of massive black holes at high 
redshift. However, our simulations are not large enough to check 
if a population of sufficiently massive black holes has formed at 
high redshift to account for observations of quasars that have been 
interpreted as implying the presence of very massive black holes 
at early times (Fan et al. 2001, 2004). Fan et al. have discovered 
quasars of magnitude M1450 ~ —27 at 2 ~ 6, albeit at a low space 
density, <^ 1.6 ± 0.5 x IQ-^^^^Mpc-^ dFan et al. 2004t . As- 
suming that these objects are radiating at the Eddington luminosity 
with an Elvis et al. d 19941 ) spectrum, and that beaming and gravita- 
tional lensing are insignificant. Fan et al. ( |2001) inferred that these 
quasars host black holes masses of ~ 1 - 3 X lO®/i"^M0. In this 
paper, we probe the mass function of 2 = 6 black holes only down 
to a space density of <^ ~ 1 x 10~*/i"^Mpc^'*, which is not sparse 
enough to compare with these data. We plan to address this problem 
by performing simulations of a larger volume than those analyzed 
here, which will probe the high mass end of the halo mass function 
in detail. 

In summary, we have presented a new model for the concur- 
rent growth of galaxies and black holes in the ACDM cosmology. 
We have previously shown that this model can successfully account 
for many observed properties of the galaxy population over a large 
range of wavelengths, from the local optical and infrared galaxy lu- 
minosity function to the number counts of submillimeter galaxies 
and the UV luminosity function of Lyman-break galaxies at red- 
shift 2 ~ 3 (Baugh et al. 2005). In this paper we have focussed 
on the properties of the black hole population that grows in uni- 
son with the spheroidal component of the galaxies. This model can 
account for a variety of observables that involve black holes, such 
as the relationship between the mass of the central black hole in 
galaxies and the properties of the bulge, the quasar luminosity func- 
tion, and the apparent "anti-hierarchical" growth of black holes. 
The model may be tested by future observations of the evolution of 
the A/bh — bulge relations and, perhaps, by the detection of grav- 
itational waves associated with the mergers of massive black holes 
that play a prominent role in our model. 
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